Preprint CAMTP/96-8 
September 1996 



Supersymmetric quantum mechanics based on higher 

excited states 

Marko Robnik[] 

Center for Applied Mathematics and Theoretical Physics, 
University of Maribor, Krekova 2, SLO-2000 Maribor, Slovenia 



Abstract. We generalize the formalism and the techniques of the supersym- 
metric (susy) quantum mechanics to the cases where the superpotential is 
generated/defined by higher excited eigenstates. The generalization is tech- 
nically almost straightforward but physically quite nontrivial since it yields 
an infinity of new classes of susy-partner potentials, whose spectra are ex- 
actly identical except for the lowest m + 1 states, if the superpotential is 
defined in terms of the (m + l)-st eigenfunction, with m = reserved for 
the ground state. It is shown that in case of the infinite 1-dim potential 
well nothing new emerges (the partner potential is still of Poschl- Teller type 
I, for all m), whilst in case of the 1-dim harmonic oscillator we get a new 
class of infinitely many partner potentials: for each m the partner potential 
is expressed as the sum of the quadratic harmonic potential plus rational 
function, defined as the derivative of the ratio of two consecutive Hermite 
polynomials. These partner potentials of course have m singularities exactly 
at the locations of the nodes of the generating (m + l)-st wavefunction. The 
susy formalism applies everywhere between the singularities. A systematic 
application of the formalism to other potentials with known spectra would 
yield an infinitely rich class of "solvable" potentials, in terms of their partner 
potentials. If the potentials are shape invariant they can be solved at least 
partially and new types of analytically obtainable spectra are expected. 
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1 Introduction 



After the classical papers of Witten (1981) and Gendenshtein (1983) the 
methods of supersymmetric (susy) (nonrelativistic) quantum mechanics have 
quickly developed and it has been realized, that (1) there exist partner po- 
tentials with precisely the same energy spectra except for the ground state 
(m = 0) (whose wavef unction <p(x) = if)o(x) is used to generate/define the 
superpotential W(x) - see below) 0, and that (2) if they are "shape invari- 
ant", their spectra and wavefunctions can be exactly and analytically solved. 
It is believed that the list of such shape invariant partner potentials is now 
complete and finite (Levai 1989, Barclay et al 1993), and therefore quite lim- 
ited in use. The research has been later further developed also in direction 
of applying the WKB methods to such classes of Hamiltonians, including 
the search for improved simple quantization conditions which would be ex- 
act in case of susy shape invariant potentials (Barclay, Khare and Sukhatme 
1993, Barclay and Maxwell 1991, Barclay 1993, Inomata, Junker and Su- 
parmi 1993, Junker 1995, Robnik and Salasnich 1996), and also in direction 
of exploring the applicability of the path integral techniques (Inomata and 
Junker 1991,1994). One of the nicest presentations of susy quantum mechan- 
ics was published by Dutt, Khare and Sukhatme (1988), henceforth referred 
to as DKS. We will use their notations. It should be mentioned at this place 
that the ideas involved behind the susy property and shape invariance were 
formulated first by Infeld and Hull (1951), where they were called the "fac- 
torization method", and these authors refer further to the related ideas in 
the works of Schrodinger (1940,1941). 



2 Generalized supersymmetric formalism 

The main point of this short Letter is to point out that the whole formal- 
ism of the susy quantum mechanics can be generalized to arbitrary higher 
excited eigenstates <j){x) = i[) m (x), m = 0,1,2,..., used to generate the 
superpotential W(x), namely 

2 the ground state energy Eq ' is missing in the partner Hamiltonian H +1 so that its 
groundstate E ( +) = e[~ } 
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W(x) 



(i) 



where (f>'(x) = d<fr/dx, \x is the mass of the particle moving in the V_ potential, 
27rh is the Planck constant and m is the quantum number equal to the number 
of nodes of the eigenf unctions ip m (x) of the starting potential V-(x). The 
energy scale is adjusted so that the (m + l)-st energy eigenvalue is exactly 
zero, Ej~^ = 0. The corresponding Hamiltonian is H = — IrJ^ + V-(x), 
and the Schrodinger equation reads 



Obviously, because 4>\x) ^ at the nodes yj, the superpotential W(x) will 
have singularities at the nodes yj, j = 1, 2, . . . , m of 0. However, this does 
not invalidate our derivation, but it merely means, as will become clear later 
on, that the partner potential generated by (ft diverges to +oo when x — > yj, 
for any j — 1, 2, . . . , m. This implies that the potential wells are well defined 
between two consecutive singularities and that they do not communicate with 
solutions in the neighbouring wells. Thus if m = we have the common case 
of usual susy potentials defined on (— oo, +oo), if m — 1 we have two sepa- 
rated potential wells, each of them on a semiinfinite domain, for m = 2 we 
have one infinite potential well on a finite domain between two nodes y\ and 
?/2, and two binding potential wells on the two seminfinite domains (— oo, yi] 
and [2/2, +oo), and so on. The (partner) potentials constructed in this way 
are nontrivial and certainly very interesting since they contribute to our list 
of solvable potentials which now becomes truly very rich and infinite in its 
contents. 

In order to make this Letter selfcontained I will build up the formalism nec- 
essary to construct the partner potentials and to define the shape invariance, 
following DKS, in order to demonstrate that the susy formalism does not 
break down anywhere on its domain of definition, and to define the language 
needed to talk about further results that I shall present in this contribution. 

First we express the starting potential V-(x) in terms of the (m + l)-st 
eigenfunction <j)(x) = if) m (x), by solving @ 



= H_<P=( 



2fi dx 2 



+ V4x))(f> 



0. 



(2) 
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which is regular everywhere, because at the nodes z/j the second derivative 
4>"(x) = d 2 (p/dx 2 also vanishes with 0. Thus the basic Hamiltonian H reads 

The two important operators are: 

and 

a n ( d V\ (a\ 



which gives 

H = A + A. (7) 
We further define the partner Hamiltonian H + and the partner potential V + 



as 



where 



h += aa+ = -tS +v+(x) < (8) 



or 



V + (x) = -V4x) + -(^-) 2 . (10) 

The potentials V + and V_ are called susy-m partner potentials. We will show 
that they have the same energy levels, except for the (m + 1) lowest states 
of V_ for which there are no corresponding states of V + , so that the ground 
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state of the latter one is Eq = E^+i- All higher states have then iden- 
tical energies. From equation fllCf) we see explicitly that at every node jjj, 
j = 1,2, ... , m of the defining eigenstate <p — 4>m the partner potential(s) 
will have a singularity of the type l/(x — yj) 2 such that V + (x) — > +00 when 
x — > yj, so that every branch of the partner potential will be confining up to 
infinity, and the solutions in various branches do not communicate. Thus for 
each m we shall find (m + 1) (branches of the) partner potentials. 

In terms of the superpotential W defined in equation (|l|) we can write 

<j>{x) = if>(-\x) = e*p(-^ f W(x)dx), (11) 

which is well defined in the definition domain of any of the branches of the 
partner potential, and obviously <ft will go to zero where W has the poles 
l/(x — yj), as it should happen. 

Some of the relationships can be rewritten/reformulated in terms of the su- 
perpotential W(x) now: 



A = -7T»T X + W{X) - 



A h d T T 7" / N 

A = 7mTx + W{x) > 



and also, the commutator of the operators A and A + is 



Further we observe 



and also 



(12) 



V±(x) = W 2 (x) ± -=W'(x), W'{x) = (13) 



„ „ 2ft dW 

v - = v - + (14) 



KM 

L ' J v^u dx v ; 
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Now we have all tools at hand to show that the susy partner potentials V- 
and V+ are isospectral except for the lowest (m + 1) states of V- which have 
no counterpart in V+, so that its ground state is Eq + ^ = 

The demonstration, following DKS, is very easy: First we find that if is 
an eigenfunction of f/_ with the eigenenergy E^~\ then Aip^ is an eigen- 
function of H + with the same energy: 



H + (A^- ] ) = AA+A^ = AH-^P = AE^^ = E^A^~\ (16) 

Now we show that this applies only to the eigenstates n higher than m, n = 
m + 1, m + 2, . . ., by considering the normalization condition, by writing the 
normalized state ip^ = C n Aip(~\ and calculating the normalizing coefficient 



ll^ +) ll = Cl < A^\A^n ] >= Cl < ^{A+A^ >= ClE^\\^-\ 

(17) 

If all ip^ > are normalized (they are certainly orthogonal, because we deal with 
one dimensional systems, where degeneracies are forbidden due to the Sturm- 
Liouville theorem (Courant and Hilbert 1968) and therefore all eigenstates 
must be orthogonal), then 

which implies that the construction succeeds only iff E^ > 0, implying that 
n > m. Thus the two Hamiltonians H_ and H + defined in (|]) and in (H) 
are isospectral except for the lowest (m + 1) eigenstates of H- which have 
no counterpart in H + . 

Counting now the eigenstates of H + from n = 0, 1, 2, . . ., where n = is 
the ground state, and n is the number of nodes of the (now also normalized) 
eigenfunction ip^\ we have 

= 4 +) = E { ~l +n . (19) 
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Of course it is easy to show that, conversely, for every eigenstate ip^ of H + 
there exists the normalized eigenstate of namely 




1 



(20) 




n = 0,1,2,... 



This completes our proof of isospectrality, generalized to the case that the 
generating function <fi of the superpotential W, defined in equation (fj), is a 
higher excited wavef unction, namely <p — V'm > wi = 0, 1, 2, . . .. As we have 
seen, the formalism of superpotential and of the partner potentials works 
everywhere except at the singularities located at the nodal points yi of </>, 
where the partner potential V+ goes to infinity as l/(x — yi) 2 , thereby defin- 
ing several branches of V + well defined on their disjoint domains of definition. 

We have demonstrated that if one of the partner systems (the Hamiltonians) 
can be solved completely (by calculating the energy levels and the eigen- 
functions), then the susy formalism enables one to solve the partner problem 
completely, following equation (Tl9|). One of the most important cases is of 
course the harmonic oscillator, which we will discuss in detail below. 

If the solutions for the two partner Hamiltonians are both unknown, then 
another approach is necessary to solve them. In case of the standard susy 
formalism with m = we have the important class of the shape invariant 
potentials. As is well known (DKS) the shape invariance of the two partner 
potentials V_ and V + is defined by 



where a is a set of parameters, a± = f(a ) and R(ai) is independent of x. 
The procedure is now (essentially embodied in the factorization method of 
Infeld and Hull (1951)) the following. Consider a series of Hamiltonians H^ s \ 
s = 0,l,2,..., where H<® = H_ and = H+, by definition 



V + (x; ao) = V-(x; di) + R(ai), 



(21) 



H (s) 



2/i dx 2 



+ V-(x; a s ) + R{ a k), 



k=l 



s 



(22) 



where 
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a B = / fl (ao)= /o. . .o/ (a ). (23) 

s 

Now compare the spectra of H^ a ' with H^ s+1 \ and find 

= -^ + V4x;a s+1 ) + Y / R(a k ), 
z ^ k=i 

# (s+1) = -^ + V+ (x;a s ) + j2R(a k ), (24) 

" k=l 

Thus it is obvious that and if are susy partner Hamiltonians, and 
they have the same spectra from the first level upwards except for the ground 
state of whose energy is 

E { s) = j2R(a k ). (25) 

k=l 

When going back from s to (s — 1) we reach = H + and = H_, 
whose ground state energy is zero and its n-th energy level being coincident 
with the ground state of the Hamiltonian H^ n \ n = 1, 2, . . .. Therefore the 
complete spectrum of if_ is 



= X>K), Et ] = 0. (26) 

k=l 



The generalization of shape invariance to the case of any m > is straight- 
forward, but it results in higher complexity and therefore it is more rarely 
satisfied by the specific systems. By repeating the above argumentation we 
reach the conclusion that, when (j2l"D is satisfied for a superpotential W with 
given m, then we cannot calculate the entire spectrum of the shape invari- 
ant potential/Hamiltonian H_, but only the subset (subsequence) of period 
m + 1 , namely 

n 

E { ~l n(m+1) = E RM, E<r) = 0, n = 1, 2, . . . . (27) 
fc=i 

In the special case m = we of course recover the formula (|26|) . For m > 
we have none example of susy-m shape invariance so far. 
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3 The example of the harmonic oscillator 

Let us consider a few examples of susy-m partner potentials, first the har- 
monic oscillator, defined by 



V-{x) = i/iaA 2 - (m + -)hcu, (28) 



shifted in energy so that 



= 0. (29) 




Introducing the natural unit of length a we can write the ground state wave- 
function as 

1 x 2 

= 7rl/4al/2 eX P(-^)> « = J" > (30) 

Defining the creation operator a + , 

" + = ^-"l + (31) 

we can write down all the eigenfunctions, in particular the (m + l)-st one, 
labelled by m, and denoted by 0, as follows 

(a + ) m 

^-(x) = ^ r ^o(x)=0(x). (32) 

Now we calculate the superpotential according to equation ([]]), by using the 
following operator when calculating <fi' = d<p/dx, obtained from equation 

^ = -(--v/2a + ), (33) 
ax a a 

and find 

W(x) = — * ( X - - . (34) 

Using the explicit solution for ip m (x), namely 
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1/2 2 

.X s , X 



*"W=^VS?J H - ( a )eXp( -2^)' (35) 

where H m (z) is the Hermite polynomial (Abramowitz and Stegun 1965), we 
obtain 



An+l(x) 1 H m+ l( 



x 



Mx) ^2(m + 1) H-m(£ 

and hence 



(36) 



W(x) = — * ( X - - . (37) 

a^2fx \a Yl m [x/oi) J 



From this and equation ([14]) we get finally the susy-m partner potential of 
the harmonic potential (p8|), namely 



K +W = i^V- (m+ ^ + ^(^±|)) (38) 

whose energy levels are the same as for the (energy shifted) harmonic oscil- 
lator, namely 

E { n +) = (n + l)hcu, n = 0,1,2,... (39) 

so that the ground state energy = hu, which is equal to the (m + l)-st 
energy level of (|2~8|). 

This result is important, because for m > it yields new interesting poten- 
tials with purely discrete spectrum, isospectral to the harmonic oscillator, 
except for the lowest (m + l)-st eigenstates. Let us look just at the few 
lowest cases. 



m = : V + {x) = -fiu 2 x 2 H — hu, 

2 2 

-r t 1 \ 1 9 2 1+ ^ 02 

m = 1 : v+vx) = -jj,uj x nuo + Too—, 

2 2 x 1 
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2^ 2 \2z 2 - 1 (2z 2 -I) 2 J' a' 

rr , s 1 2 2 5 + /, 6 12 \ 

m = 3: V+(x) = ~nux 2 - -Tiu + 1 + — - + 



2^ 2 V 2z 2 -3 (2z 2 -3) 2 y ' 

The case m = is the usual susy-0 case showing just that the 1-dim har- 
monic oscillator potential is indeed susy-0 shape invariant. But m = 1 gives 
a new example, which nevertheless is well known as the radial problem of 
the 3-dim harmonic potential, which is thus susy-1 partner potential of the 
1-dim harmonic oscillator potential. It has only one singularity at x — 0. 
Next, for m = 2 in the above list, we see the first new nontrivial example, 
of a specific rational potential which is susy-2 partner potential of the 1-dim 
harmonic potential. It has singularities at the two nodes x — y\ — —a/\r2 
and x = y2 = +a/y/2, of the type l/(x— Uj) 2 - Therefore it has three branches 
(ranges), namely (— oo,yi], [?/i,?/ 2 ], and \y 2 ,+oo). The spectrum is identical 
in each of them. Further, in case m = 3, we have three singularities at the 



nodes where x is equal to y\ = —aJ3/2, y 2 = 0, and y 3 = —y\ = +ay3/2, 
and thus we have two independent different potentials within the two ranges 
[7/2, 2/3] and [2/3, +00). (The other two ranges confine the potentials which are 
equivalent due to the evenness of V + (x).) 

For higher m we get new classes of interesting rational potentials, all of 
them being isospectral to the harmonic oscillator except for the lowest m + 1 
eigenstates of the latter which are missing in the partner potentials. For each 
m we have rational potentials with m + 1 branches, defined by the m nodes 
r/j, j = l,...,m. Since the Hermite polynomials B. m (z) are even or odd 
functions of z, depending on whether m is even or odd, the superpotential 



(37) is always odd function of x and therefore the susy-m partner potential 



V + (x) in equation (BSf) is always even function of x with m singularities. 



Asymptotically when |x| — ► 00 the potential behaves still as harmonic quadratic 
potential with the leading term ^fiu 2 x 2 , which is true for any m, as can be 
shown using the asymptotic properties of the Hermite polynomials, as de- 
scribed below. 

The limiting (semiclassical) behaviour of the potential V + when h — > is 
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interesting. It implies that z = x/a = x^Jfiu/h tends to +00 and therefore 
from the asymptotic properties of the Hermite polynomials H m+ i (z) /K m (z) —. 
2z, for z — ► +00, we conclude 

1 11 

V+(x) — > -[iuj 2 x 2 — (m — — > -/ic<j 2 £ 2 when % — > 0. (41) 

Thus the semiclassical limiting form of all these potentials is just the har- 
monic quadratic potential, meaning that all the rational potentials in (R^) 
all have zero classical limit. From (|38|) it is clear that the harmonic oscilla- 
tor potential is not susy-m shape invariant, except for m = 0, which is the 
familiar case of shape invariance (see DKS). 

4 Discussion and conclusions 

Using the same formalism applied to known solvable potentials for various m 
we can systematically construct the vast class of new potentials which will 
be isospectral to each of the known solvable potentials, almost all of them 
being susy-0 shape invariant, and listed in DKS. 

Finally we can state the result which can be easily verified in a straight- 
forward manner (we omit the derivation due to the lack of space here), that 
the Poschl- Teller type I potential is the susy-m partner potential of the infi- 
nite potentiall well for any value of m > 0. At present we do not know any 
specific cases of susy-m shape invariance with m > 0, and also have no fur- 
ther calculations for susy-m partner and solvable potentials, which remains 
as the future project. 
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